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Abstract 

In this paper we consider mixed two-loop electroweak corrections to the top 
quark propagator in the Standard Model. In particular, we compute the on-shell 
renormalization constant for the mass and wave function, which constitute building 
blocks for many physical processes. The results are expressed in terms of master 
integrals. For the latter practical approximations are derived. In the case of the 
mass renormalization constant we find agreement with the results in the literature. 


PACS numbers: 14.65.Ha 12.38.Bx 


1 Introduction 


The outstanding precision reached at the CERN LEP and SLAG SLC triggered many 
higher order calculations in the Standard Model (SM) of particle physics. In particular, 
it happened for the first time that the experimental results were sensitive to the weak 
part of the SM. Since at LEP and SLC real top quarks could not be produced, the main 
emphasis of the theoretical investigations was put on processes with light quarks as exter¬ 
nal particles. Heavy particles like the top quark only appeared virtually in intermediate 
states. In many applications the masses of the light quarks can be neglected as com¬ 
pared to the other mass scales, which results in a signihcant simplihcation of the resulting 
mathematical expressions. 

In a future International Linear Collider (ILC) [1] the center-of-mass energy is high 
enough to produce top quarks. The expected experimental precision requires on the 
theoretical side the inclusion of higher order corrections — both for the QCD and the 
electroweak sector of the SM. This is particularly true for the threshold production of top 
quark pairs, where the theoretical uncertainties of the second order QCD corrections are 
still signihcant [2]. Thus, next to third-order QCD calculations (see, e.g.. Refs. [3-6]), 
also electroweak corrections have to be investigated. 

In this paper we focus on the two-loop mixed electroweak/QCD corrections to the 
on-shell top quark propagator. Some sample diagrams are shown in Fig. ^ The top 
quark mass renormalization constant has been evaluated in Ref. [7]. The contribution 
from the scalar bosons has been considered in [8,9]. In this paper we conhrm the results 
of Ref. [7]. Furthermore, we compute the wave function renormalization constant to 
order aag- Practical approximations are derived for the diagrams with internal W and Z 
bosons^ and both the exact expression and handy approximations are computed for the 
Higgs boson mass dependence. 

The inverse fermion propagator can be decomposed as 

S~^{q) = Z^{j - Z^m + i Si(g) -F mSs) L + Z^ - Z^m + i -F mSs) R , (1) 

with R = (1 -|- 75)/2 and L = (1 — 75 )/ 2 . m is the quark mass and Z^ and Z^’^ are the 
mass and (left/right) wave function renormalization constants, respectively. The functions 
Si, Si; and S 5 in Eq. (|T]) result from a convenient decomposition of the self energy S(g) 
given by 


S(g) = ^(RSi(gQ + LSi(gQ)+mSs(gQ. (2) 

In the on-shell scheme one requires that S~^{q) vanishes for q^ = m? which leads to 
the following condition for the on-shell mass renormalization constant 

27 = 1 + + i (SL(m^) + . (3) 

^Here and in the following we consider the W and Z boson always in combination with the corre¬ 
sponding Goldstone bosons. 
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In these formulae it is understood that only the real part of the self energy functions is 
taken. 

In the practical calculation it is convenient to apply projectors in order to arrive at 
scalar momentum integrals. This is achieved via 



Figure 1: Sample diagrams contributing to the top quark propagator up to 
order aag- Next to the Higgs boson also the gauge bosons W and Z and the 
corresponding Goldstone bosons can be exchanged. 
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In these formulae it is understood that the lower-order results are expressed in terms of 
the bare parameters. In what follows we compute and for the top quark 

neglecting the masses of the light quarks. Furthermore, we take the Cabibbo-Kobayashi- 
Maskawa matrix to be diagonal. 

The outline of the paper is as follows: In the next Section we consider the classes 
of Feynman diagrams which are relevant for our calculation and discuss the reduction 
to a basic set of master integrals. The results for the master integrals appearing in our 
calculation are discussed in Appendix O In Section El we present the renormalization 
constant for the top quark on-shell mass. In particular we consider the relation between 
the MS and pole mass and compare our results with the literature. In Section 0] we move 
on to the wave function renormalization constants. Finally, our hndings are summarized 
in Sectional which also contains the conclusions. 


2 Reduction to master integrals 

At order aag one has to consider about 35 Feynman diagrams contributing to the 
fermion propagator (cf. Fig. ^ for some sample diagrams). After the application of the 
projectors the external momentum is set on the mass shell of the heavy quark, which leads 
to integrals containing two scales, the quark mass and the boson mass. In this Section 
they are denoted by m and M, respectively. 

We generate all one-particle irreducible Feynman diagrams contributing to the fermion 
propagator with QGRAF [10]. The application of q2e and exp [11,12] identihes the topol¬ 
ogy of the individual diagrams, adopts the notation and transforms the expressions into 
FORM [13] notation. 

It is convenient to map the QGRAF output for each diagram with the help of q2e and 
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Figure 2: Graphical representation of the integral classes defined in Eq. fTTjl. 
Solid and dashed lines carry mass m and M, respectively. Curly lines are 
massless. 


exp to one of the following four classes of integrals^ 


B'^(ni,n2,n3,n4,n5) = 
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(7) 


where d = 4 — 2e is the space-time dimension. The corresponding graphical representation 
can be found in Fig. |21 The integrals are defined in Minkowskian space and the ie 
prescription is understood. 

Within FORM we apply the projectors, identify the external momentum with the top 
quark mass and decompose the numerators in terms of the denominators. This leads 
to a large number of scalar integrals which differ from each other by the power of the 
individual propagators. 

A conventional way to reduce an arbitrary integral of a certain kind to a small set of 
so-called master integrals is based on integration-by-parts [14], which provides relations 
between several integrals of different complexity. The proper combination of these rela¬ 
tions leads to new ones, so that the iteration of this procedure can be used for a systematic 
reduction of an arbitrary integral to a small set of master integrals. 

^For the application at hand it is possible to work with a smaller set of integral types. E.g., in the 
case of we have 715 < 0 for the self energies considered in this paper. However, in view of a future 
application to on-shell vertices it is advantageous to consider a more general set-up. 
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In Ref. [15] an algorithm has been formulated that performs automatically the afore¬ 
mentioned reduction for a given set of recurrence relations. Currently several implemen¬ 
tations of this algorithm exist. However, to our knowledge, only one program is publicly 
available, AIR [16]. AIR is written in MAPLE, which certainly constitutes a serious restric¬ 
tion for large-scale problems like, e.g., four-loop vacuum integrals [17-20]. Nevertheless, 
for the problem at hand AIR is well suited to perform the reduction. It is straightforward 
to compose an interface which, for a given class of Feynman diagrams (cf. Eq. ((7j)), pro¬ 
duces the corresponding integration-by-parts relations, passes them to AIR and transforms 
the output into a table that can be read into FORM. These tables can be used to express 
each integral occuring in our expression in terms of a few master integrals. 

At this point a comment concerning AIR is in order. Our experience with AIR shows 
that there can be situations where the set of master integrals is not minimal, although the 
complete set of recurrence relations has been provided. Consider, e.g., the integrals 
with M = 0. It is well-known that only three master are needed for the computation of 
the two-loop pure QCD corrections to However, the naive application of AIR leads 
to four master integrals. A straightforward inspection of the involved integrals makes it 
possible to relate the additional master to the known ones. The same is also true for our 
types of integrals. 

For the diagrams where a neutral boson is exchanged one requires altogether nine 
master integrals. Six of them either contain only one mass scale, are vacuum diagrams, 
or consist of a product of two one-loop integrals. They read 

Hi = H+(l, 1, 0, 0, 0), H2 = H+(0, 0,1,1,1), Ha = H+(l, 1, 0,1, 0), 

Yi = Y^(l, 1,1, 0, 0), Y 2 = Y^(l, 0, 0, 0,1), Ya = Y^(l, 1, 0, 0,1). (8) 

Explicit analytic results are given in Appendix O The remaining three master integrals 

H 4 = H+(l, 0, 0,1,1), Hs = H+(2, 0, 0,1,1), Y 4 = Y^(l, 1,1, 0,1), (9) 

are less trivial. Analytic expressions for H 4 , H 5 and Y 4 can be found in Ref. [7]. More 
precisely one has 


H 4 Joi 2 (l, 1, 1, Af^) , 

H5 <-»• ^012(1, 2 , 1 , m^, M^), 

T4 ^ 1, 1, 1) , (10) 

where the integrals Joi2(l, 1,1, ''^^5 and Joi2(l, 2,1, m^, Af^) are given in Eq. (3.20) 
and I 5 1,1) In Eq- (3.26) of Ref. [7]. We have checked all master integrals by 

considering their evaluation in an asymptotic expansion around the three kinematical 
regions m -C Af, m ^ M, and m 3> Af. Altogether we computed up to 16 expansion 
terms and found complete agreement with the results in the literature. As can be deduced 
from the results of Appendix [Cl where more details are provided, the inclusion of about 
hve expansion terms in each region provides jointly a good approximation over almost the 
whole range in m/M. 
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In the case of charged boson exchange one gets in addition two simple master integrals 
11^1 = 11^J(1,1,1, 0, 0), W2 = 11"J(1,1, 0, 0,1), (11) 

and hve two-scale integrals 

H, = H+{0, 1,1,1, 0), Hr = H+{0, 1,1,1, -1), Hs = H+{1, 1, 0,1, 0), 

hh3 = irj(0, -1,1,1,1), W4 = lhj(l, 1,1, 0,1). (12) 

As we will see in Sections 0] and lU for the physical applications of this paper an expansion 
for m 3> M of the integrals in Eq. m is sufficient to obtain hnal results which in the 
physical region are equivalent to the exact expressions. 


3 On-shell mass renormalization constant 


In this Section we discuss the results for the on-shell mass renormalization counterterm. 
The QCD corrections up to three loops can be found in Refs. [21-25]. The one-loop 
electroweak and two-loop mixed corrections for light quarks can be found in Refs. [26] 
and [27], respectively. In this case it is sufficient to evaluate the limiting behaviour for 
rrig -C (where M represents a boson mass). The corrections of order aas for the 
top quark have been considered in Ref. [7-9]. In a recent paper [28] the two-loop relation 
between a minimal subtracted and the on-shell mass has been considered in a more general 
framework. However, the masses of the vector bosons have been neglected. 

The relation between the bare mass, m°, and the one dehned in the MS and on-shell 
scheme, mt and rrit, is given by 

= Z“m,. (13) 


In order to discuss the result for the mass renormalization constant it is convenient to 
consider the hnite ratio 
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with an analog separation for and We furthermore introduce the notation 


Vh 

Lh 


rrit _ _ rrit 

“ ~ ’ 

= — ln(?/^) , Lyy = — ln(j/^) , Lz 


= - HVz) > 



The renormalization constant in the MS scheme is given by (see, e.g., Ref. [7]) 
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with at = l/{2swcw), Vt = (1/2 - 4s^/3)/(2svvCw) and cw = a/1 - = Mw/Mz- For 

convenience, in Eq. (ED the contribution from the tadpole diagrams is displayed separately 
in the second {0{a)) and last line {0{aas)). Furthermore, all term proportional to W 
originate from tadpole diagrams. 

At this point a comment concerning the various gauge parameters is in order. In the 
electroweak sector we adopt Feynman gauge for the W and Z boson, however, we allow 
for a general gauge parameter ^ for QCD defined via the gluon propagator 


^r(^) 






—z ■ 




+ ie 


( 18 ) 


On general grounds the on-shell mass and also Zm has to be independent of ^ which serves 
as a welcome check for our calculation. 

In Eq. ED the two-loop expression for z^ is split into contributions induced by the 
Higgs, W and Z boson (including the corresponding Goldstone parts), the photon (A) 
and the tadpole diagrams. In the following we present analytical results for the individual 
contributions. As we will see in the discussion below it is sufficient to consider the limit 
yw —^ C )0 and yz —>■ oo in order to obtain agreement with the exact result below the per¬ 
cent level. For this reason we show only the corresponding analytical expressions. Since 
























the Higgs boson mass is still unknown we present the exact result, but also handy expan¬ 
sions in the three limits uh ^ Uh ^ ^ and yn —^ C) 0 . Furthermore, for completeness 
also the tadpole result is listed, which can be extracted from Refs. [8,27,29]. 

Compact expressions for the Higgs boson contribution to and z^™^^ expressed 

in terms of (known) master integrals are given in Eqs. (jTfj) and pHjl in Appendix fAl 
where also the results for the //-dependent terms and can be found. 

The expansions in the physically interesting regions read for the one-loop results 
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with l — The subscripts 0, 1 and oo indicate the cases rrit Mh-, rrit ~ Mh-, 

and Wit 3> Mh,w,z, respectively. The coefficients in front of are given by 
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The corresponding two-loop expressions are slightly more lengthy but still rather com¬ 
pact. They read 
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and C, is Riemann’s zeta function. The dehnition of Lsi(z) is given in Eq. (14211 . The /i 
dependence is determined through 
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(a) 


(b) 


Figure 3: (a) One-loop and (b) two-loop corrections to Zm as a function of 
^/Vh — Mn/mt- The solid (coloured) lines correspond to the highest available 
order for each case. The dotted curves show lower-order results and nicely 
demonstrate the convergence. The exact result, which is plotted (in black) 
over the whole l/y^ range, can be distinguished from the approximations 
only in small gap around l/y^ ~ For the renormalization scale /i^ = 
has been chosen. In the plots the contributions from the tadpole diagrams are 
not included. 


Note that the result for in the above expressions can easily be extracted from the 
two-loop pure QCD result [21]. 

Let us next discuss the numerical consequences of our result and compare the exact 
expressions with the compact expansions. Actually the exact result from the W and 
Z boson contribution is very well reproduced both by the large-top mass limit or the 
expansion around rut ~ Mw/z- At the physical values for the particle masses we hnd a 
deviation from the exact result for the W and Z contribution below the percent level, 
which justihes the use of the expansion for the numerical evaluations. 

The result for the one-loop coefficient of a/(7rs^) is shown in Fig. EKa) for /i = rut as 
a function of l/yn = Mn/mt where the tadpole contributions [9], which are numerically 
quite large, have been subtracted for Feynman gauge. Furthermore the following input 
values have been chosen 

TTit = 165 GeV, M\y = 80.425 GeV, Mz = 91.19 GeV, c\y = M]yIMz- (25) 

Next to the (black) solid line which includes the exact result for the Higgs mass dependence 
and the large-mt results for the W and Z contributions we also show the expansion terms 
in the three kinematical regions (cf. Eq. ()19j) i as solid lines. The lower-order results are 
plotted as dotted lines in order to demonstrate the convergence of the approximations. It 
can be seen that over almost the whole range of y^j the expansion terms provide a very 
good approximation to the exact result, except for a small region with y^ zz 0.5... 0.7 
which corresponds to Mh ^ 250... 300 GeV. We want to mention that in Fig. Eta) also 
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the result of Ref. [7] is shown which in contrast to ours also takes into account the exact 
dependence on and Mz- No visible effect is observed. 

In Fig. Efb) we show the two-loop coefficient of aas/^n'^s'^) as a function of 1/Fir, 
where again the tadpole contribution of Eq. m is subtracted. The result containing 
the exact Mh dependence (cf. Eq. 0121)) is plotted together with the expansion results 
of Eq. ()21|1 . In addition to the highest expansion terms we show as dotted lines also the 
lower-order ones. One can see how the approximations nicely improve while including 
step-by-step the higher order terms. Note that at two-loop order it is not possible to 
separate the result given in Ref. [7] into the contributions from the individual bosons and 
the tadpole contribution. Thus in order to check our results we again include the exact 
results from Ref. [7] in our plots (after subtracting the tadpole contributions). Again no 
difference is visible. In fact, adopting from Eq. m the large-m* terms we obtain at the 
physical values of Mw and Mz agreement with the exact result below the per cent level. 
Equivalently we can use the expansion around frit ~ M]y/z and get results with the same 
level of accuracy. 

The excellent description of the exact result by the expansion terms together with 
the relative simplicity of the results in Eq. m provides sufficient motivation to apply 
in the next section the same approach to the wave function renormalization constant. In 
particular this means that both exact results and expansions are considered for the Higgs 
boson contributions and approximate formulae for large top quark mass are derived for 
the remaining parts. 


4 On-shell wave function renormalization constants 


In this Section we consider the wave function renormalization constants for the top 
quark dehned through Eq. 0. Since and contain infra-red divergences which 

only cancel when considering a physical quantity, it is not possible to form a hnite ratio 
analog to in Eq. (uni)- Thus, in the following we consider the divergent contributions 
and the hnite parts separately. Furthermore, we switch to the vector and axial-vector 
contribution using the formulae 


rzV,OS 

— 

1 

2 

(zf + Z^’°® 

ryA.O?, 

— 

1 

2 

^^R,OS _ ^L,OS 


(26) 


where a non-zero contribution to only arises from the W- and Z-boson. 

In contrast to Z°® the wave function renormalization needs not to be gauge parameter 
independent since it does not pose a physical quantity. As for the mass renormalization 
constant, we choose Feynman gauge for the electroweak part but allow for arbitrary ^ in 
the QCD sector. We observe that ^ drops out in the case of Z^’°® up to the two-loop 
order, which is in analogy to QCD where only the three-loop result starts to depend on 
^ [30]. On the other hand is ^ dependent starting from order aag- 
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In order to present the results in a compact form we introduce the notation 


ryX,OS 

^2 


a 


1 + + 


TT 


2,X 


a 

TT 


few 

J2,X 


+ z'. 


2,X 


+ 


aa. 


Cf 


TT^ 


mix Lmix 

92,X , '^2,X , ^mix 
19 I : ^2,x 


(27) 


with X = V, A and the analog splitting as in Eq. m- In particular, the results for the 
individual coefficients are split according to the dependence on the renormalization scale 
and decomposed into contributions originating from the Higgs, W and Z bosons and the 
photon^ where for W and Z again the large-m* limit is adopted. 

In contrast to the quotient Zm, where it was found convenient to use the MS-mass as 
expansion parameter, for Z 2 ’ the natural choice is to express the results in terms of the 
on-shell mass. Thus, in analogy to Eq. (CHD we introduce the following notation: 
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(28) 


For the numerical value of the top quark mass we use = 175 GeV. 

In the main text we again list the expansion terms and relegate the exact expression 
for the Higgs boson contribution to the Appendix (cf. Eq. (I4f)|l b For the QCD result 
and the pole part of the one-loop electro-weak corrections we obtain the following exact 
expressions (see, e.g.. Ref. [31]) 
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The hnite one-loop contributions read 




fA,ew _ 

h,v 

fW,ew _ yw 


/ Z.ew ^ 'z 

2 A 


1 

16 ’ 


(29) 


^The tadpole diagrams do not contribute to the wave function renormalization. 
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and the scale-dependent terms are given by 
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The donble-pole parts of the two-loop resnlt are still qnite compact and can be cast 
into the form 
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whereas for the single poles we obtain the expansions^ 
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'‘Note, that the exact expressions for the Higgs boson results are given in Eq. dSU- 
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The scale dependence is ruled by the following coefficients: 
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Let us in the following discuss the results for the finite contribution to The result 
from the Higgs boson exchange expressed in terms of master integrals is given in Eq. (H) 
of Appendix El The expansion terms read 
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The results covering the scale dependence are given by 
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Figure 4: (a) One-loop, (b) 1/e pole and (c) constant part of the two-loop 
corrections to as a function of l/yn = Mn/mt- The solid (coloured) 

lines correspond to the highest available order for each case. The dotted curves 
show lower-order results and nicely demonstrate the convergence. The exact 
result, which is plotted (in black) over the whole l/yn range is only visible in 
a small gap around l/yn ~ 2. 
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In Fig. m we discuss our analytic results in numerical form. Fig. EJa) shows the hnite 
part of the one-loop contribution where the exact result is represented by the solid line 
and the expansions are plotted as dotted curves. Similarly to Zm (cf. Fig. El) it can be 
seen that the dotted lines nicely converge to the exact curve after including successively 
higher order expansion terms. As one can see, after taking into account the result from the 
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regions ^ 0,1 and cx) there remains only a quite small range for yu {^/vh ~ 1.5...2.0) 
where the (black) solid curve is still visible and the simple expansions fail to provide good 
approximations. The situation is very similar for the divergent 1/e and hnite two-loop 
contribution which are shown in Figs. Hb) andlUc), respectively. Thus we can conclude 
that it is possible to avoid the use of the quite complicated exact expressions for 
but to adopt an adequate simple expansion. For Higgs boson masses outside the range 
Mh ~ 250 ... 300 GeV one can simply use the corresponding formula from Eq. dSil) while 
for 250 GeV< Mh < 300 GeV a straightforward interpolation provides a sufficiently good 
approximation. 

It is interesting to mention that the result obtained in the gaugeless limit approxi¬ 
mates the full result within approximately 20% accuracy for Higgs boson masses between 
100 GeV and 800 GeV. For Zm the situation is similar once the tadpole contributions are 
discarded. In case the latter are included the relative deviation between the full result 
and the one in the gaugeless limit becomes smaller [9]. 


5 Conclusions 

The renormalization constants constitute building blocks for the evaluation of quantum 
corrections to various processes. In this paper the on-shell top quark mass and wave 
function renormalization constants have been considered in the SM up to order aag. The 
inclusion of electroweak effects introduces a further scale into the problem, as compared 
to the QGD or QED corrections, which makes the calculation of the integrals signihcantly 
more complicated. 

We expressed and as a linear combination of a handful master integrals 

which are known analytically, however, contain quite involved functions. For the compli¬ 
cated two-scale master integrals we applied the powerful method of asymptotic expansion 
in three different kinematical regions, which leads to power expansions multiplied by sim¬ 
ple logarithms. We could reproduce the result for available in the literature. The 
expression for is new and constitutes a building block, e.g., in the mixed elec- 

troweak/QGD corrections for top quark pair production at threshold. We checked that 
both for and the simple expansions agree very well with the exact result. 

In particular, it has been shown that the expansion for large top quark mass leads to 
compact formulae which approximate the exact results quite precisely — almost up to 
the point rrit = M with M = Mw, Mz of Mh- As far as the Higgs boson mass dependence 
is concerned, also the expansions around rrit = Mh and m* Mh have been considered. 
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A Analytic results for and 

The exact dependence on yn of the one-loop electroweak contribution reads 
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where the function Bp{ni,n 2 ) corresponds to the one-loop on-shell integral with a internal 
top quark and Higgs boson line 
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where we have included the order e terms which are needed for the two-loop expressions, 
and Tj is a shorthand notation for {y = m/M)\ 
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At two-loop order it is convenient to write the finite quantity form 

n,mix _ n,mix , J_n,mix 

where corresponds to the generic two-loop result and origins from countert¬ 

erm contributions and products of one-loop results. We obtain 
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with y^r, Vh Lfj^t are defined in Eq. (ITK|l . The poles which are still present in 
and cancel in the proper sum. 

Finally, the one- and two-loop coefficients determining the fi dependence are given by 
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B Analytic results for and 

The one-loop contribution from the Higgs boson exchange diagram to Z^y'^ reads 
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which — after inserting the results for 5,(ni,n2) from Eq. ()39j) . with m = rrit and 
M = Mh, and expanding in e — can be cast in the form 
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1 ^ ^ ^ I 

where yn and Lh are dehned in Eq. (I2HD and Kh = —— -^ 


At two-loop order we again split into two parts. 
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where corresponds to the generic two-loop result and 5^’™'^ origins from counter¬ 

term contributions and squared one-loop results. We obtain 
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(50) 


Note that the individual contributions ^4^’™“ and still depend on the QCD 

gauge parameter ^ which cancels in the proper sum. For simplicity these terms have 
already been omitted in Eqs. (gni) and dsn. 
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From the formulae ()49|1 and CT it is straightforward to extract exact expressions for 
the pole parts which are given by 


92,V 


7 

'h,v 


3y 


w 


64 


Vw 


In {Kh 


yfi 


1-Ay 


vl 


32 


Lh + 


In {Kh 


H 


11 


7 


+ 


128(1 - Ayl) 


9L 


H 


128 64(1 - Ayjj) 

7 


9Lh 

+ + 


256yj, 

9 


6Ayjj 128?/ 


H 


16^1 - Ay 


128 


(51) 


where is dehned after Eq. (HZD- 

For completeness we want to provide the exact result covering the renormalization 
scale dependence for The corresponding coefficients are dehned in analogy to the 

ones for Zm in Eq. (USD and read 
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C Master integrals 

In this Section we discuss the master integrals occuring in our calculation. The master 
integrals needed for the neutral boson exchange which reduce to products of one-loop 
integrals, to two-loop one-scale integrals or to two-loop vacuum integrals are given by 
(see, e.g.. Refs. [32,33]) 
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(53) 


with y = m/M and 


= d'i - 7^W4 - ^Lsi (ti) , (54) 

2y V r 

where 4/i, ti and ^2 are dehned in Eqs. dUl) and (EH)- 

As already mentioned in the main text, analytical results for the complicated two- 
scale master integrals can be found in Ref. [7]. We refrain from repeating them here. 
Instead we perform expansions in the limits ?/ —> 0, ?/ —> 1 and y —>■ 00 which we derived 
independently with the help of asymptotic expansions [34]. In many phenomenological 
applications it is advantageous to use the handy approximation formulae in favour of the 
complicated exact expressions. 

The expansions in the three different limits require a strategy of its own. In the case 
of a large boson mass M one obtains, next to on-shell, also vacuum integrals up to two 
loops (both for charged and neutral boson exchange) which are, e.g., implemented in 
MATAD [35]. The subdiagrams contributing in this limit can be obtained in a completely 
automated way with the help of exp [11,12]. 

In the limit y ^ 1, the case of a neutral boson exchange reduces to a simple Taylor 
expansion. The resulting integrals are well studied within QCD and documented in a 
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program code [32], The expansion for y ^ 1 oi the diagrams involving a W boson is more 
complicated, which is due to the appearance of additional massless particles, the bottom 
quarks. We did not perform the calculation of the diagrams in this limit since for the 
practical applications only the limit mt 3> Mw is needed. 

Finally for the case y —> cx), which is the limit of main interest, a careful inspection of 
the regions [36,37] contributing to the integral under consideration is in order. 

In the following we present a pedagogical example which illustrates the procedure in 
more detail. Let us consider the master integral W dehned through 


W = 


o'ie'jE 




(^^d/ 2)2 j (^2 + 2kq) {P - 2lq) {{k - If + 2q{k - /)) {P - M^) 


(55) 


In the limit y —0, i.e. m M one has to consider the cases (i) \k\ ~ |/| ~ M, (ii) 
\k\ ~ m, |/| ~ M, (hi) \k\ ~ |/| ~ m, and (iv) \k\ ~ |/| ~ M but with \k — l\ ~ m, which 
we denote as hard-hard (HH), soft-hard (SH), soft-soft (SS) and hard-soft (HS) region. 
In each region one has to perform a simple Taylor expansion of each propagator in the 
small quantities. E.g., if k is hard the propagator 1/(JP + 2kq) is expanded in 2kq/k‘^ 
since we have q^ = mP. After adding the contributions from all regions one hnally obtains 
the result for W 


W.O = E4(HH) + F4(SH) + ^(SS) + F4(HS) , 
where the integrals in the individual region have the form 

f {-2kqY^ {-2lqY^ 


(56) 


F 4 (HH) 

n(SH/HS) 

F 4 (SS) 


V c™ 

ni,...,n 5 


E 

ni,...,n 5 


^SH/HS 

Wi...n5 


C: 


SS 


'y ^ ''ni...n5 

ni,...,n 5 


(A;2)^3 ((A; - /)2 _ M2)^4 (/2)n5 ’ 

dWl {-2kqY^ {-2klY^ {-2lqY^ 

(A;2 - M2)^4 (/2 + 2lqY^ ’ 

dWl (A;2)”i (/2)^2 

(A;2 + 2kqY^ {P - 2lq)^^{{k - If + 2q{k - l)Y^ 


.(57) 


We have used partial fractioning and rearranged the different propagators in a conve¬ 
nient way. The indexes rii to run over appropriate integer values, and the C\ 


ab 

ni...n5 


are 


shorthand notations for a set of functions depending on both scales M and m and the 
dimension d. The hrst four terms of 14 expanded for y ^ 0 read 
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(58) 
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where we refrain from listing the contributions of order e. 

In the region y cx), where M m a straightforward inspection reveals only two 
different ranges for the integration momenta k and 1. In the hrst one (hh) both momenta 
are of order m, which now is hard; for the second one (hs) k is hard and I is soft, i.e. of 
order M. Thus can be written as 


n = l4(hh) + F4(hs), 


(59) 


where the following integrals occur 
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+ 2/cg)”3 [P — 

(60) 


Note that the collection of integrals in y 4 (hh) and the ones in 14 ( 88 ) are equivalent, and, 
at the same time, equal to the ones one would need to solve in order to get an arbitrary 
two loop 0{ag) strong QCD correction. The hrst seven terms of I 4 expanded for y —> cx) 
read 
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(61) 


Finally, let us consider limit the y ^ 1. Here the asymptotic expansion leads to only 
one region which corresponds to a naive Taylor expansion in the quantity A = {mP — 

This leads to the integrals of the type 




EC 


_ d^^kdH (A)^i-^ _ 

(pp + 2kq) {P — 2lq) {{k — ly + 2q{k — /)) {P — m?Y^ 


(62) 


The hrst six terms in the expansion are given by 
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\m^ J 

2e2 ^ 


>^ 4.1 = 





(63) 


with yi = 1 — l/y'^ = ^/rn? and S 2 as defined in Eq. (1^ . 

We proceeded in an analog way for the remaining two integrals in Eq. Fnrther- 
more, for the seven integrals in Eqs. dm and m we evalnated the phenomenological 
limit y ^ 00 . However, we refrain from listing the results explicitly. 
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